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We present a mathematical formalism for the description of unrestricted quantum walks with 
entangled coins and one walker. The numerical behaviour of such walks is examined when using a 
Bell state as the initial coin state, two different coin operators, two different shift operators, and one 
walker. We compare and contrast the performance of these quantum walks with that of a classical 
random walk consisting of one walker and two maximally correlated coins as well as quantum walks 
with coins sharing different degrees of entanglement. 

We illustrate that the behaviour of our walk with entangled coins can be very different in compar- 
ison to the usual quantum walk with a single coin. We also demonstrate that simply by changing 
the shift operator, we can generate widely different distributions. We also compare the behaviour 
of quantum walks with maximally entangled coins with that of quantum walks with non-entangled 
coins. Finally, we show that the use of different shift operators on 2 and 3 qubit coins leads to 
different position probability distributions in f and 2 dimensional graphs. 
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I. INTRODUCTION 

In recent years interest in the field of Quantum Walks 
has grown hugely, motivated by the importance of Clas- 
sical Random Walks in Computer Science, as well as the 
advantages that Quantum Walks may provide us with 
when compared to their classical counterparts. 

Classical random walks are a fundamental tool in Com- 
puter Science due to their use in the development of 
stochastic algorithms pj. In both theoretical and ap- 
plied Computer Science, stochastic algorithms may out- 
perform any deterministic algorithm built to solve cer- 
tain problems. A notable example is that of the best 
algorithm known so far for the solution of 3-SAT, a fun- 
damental problem in Computer Science which relies on 
random walks techniques 0- 

So random walks are important elements of Com- 
puter Science. Additionally, the recent development of 
Quantum Computation and Quantum Information has 
revealed that the exploitation of inherently quantum me- 
chanical systems for computational purposes leads to a 
number of significant advantages over purely classical sys- 
tems. Thus it is reasonable to expect that the study of 
random walks using quantum mechanical systems may 
prove fruitful. The basic properties of two kinds of quan- 
tum walks have already received a good deal of attention: 
continuous (0-0) an d discrete (01 " E3) quantum 
walks. In this paper we shall focus on discrete quantum 
walks. 

Quantum walks are expected to play a major role in 
the field of Quantum Algorithms. A number of ben- 
efits of employing such walks are already known. In 
|20|. Childs et al have built a quantum algorithm based 
on a continuous quantum walk to solve the problem of 
traversing a graph G with an exponential number of ver- 
tices under a vertex connectivity arrangement. It was 
proved in pol | that traversing such a graph with their 



quantum algorithm is exponentially faster than doing so 
with any classical algorithm. Additionally, it has been 
proved by Ambainis et al that quantum walks on a line 
spread quadraticall y fa ster than classical random walks 
|17j . while Kempe [lj| proved that the hitting time of 
a quantum walk on the hypercube is polynomial in the 
number of steps of the walk (the classical counterpart 
takes exponential time). Also, Shenvi et al showed 
that quantum walks can be used to produce a quantum 
search algorithm. Finally, an algorithm based on quan- 
tum walks to solve the problem of element distinctness 
can be found in plj . An excellent summary of the ba- 
sics of quantum walks can be found in |22j | , and a com- 
pendium of algorithmic applications of quantum walks is 
presented in |23 |. 

A discrete quantum walk is composed of two physical 
systems: a walker and a coin (a detailed explanation of 
these two systems is provided in Section III). The prop- 
erties of quantum walks applying multiple quantum coin 
operators C[24| - [26J) as well as decoherent coins ([23] - 
[2!j) on a single walker have been extensively studied. 

However, the use of entanglement in quantum walks 
is less well explored. A discussion on discrete quantum 
walks using non-separable evolution operators and its ef- 
fects on the stardard deviation of resultin g pr obability 
distributions is given in [3^. followed by |25j ] where a 
more exhaustive study on non-separable operators is pro- 
vided. In |3flj . Du et al proposed an implementation of 
a continuous quantum walk on a circle, and numerically 
showed that entanglement in the position states shapes 
the position probability distribution. More recently, a 
discussion concerning models of a quantum walk on a 
line with two entangled particles as walkers is provided 
in |34| . A study of entanglement between coin and walker 
in quantum walks on graphs is given in |3lj| . along with 
a g eneralization of the quantum walk algorithm from 
|20j |. In (3^ the authors analyse the relation between 
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coin entanglement and the mean position of the quan- 
tum walker for 3 and 4 qubit coins. Finally, in 33], Abal 
et al have quantified the entanglement between walker 
and coin generated by the shift operator in a single coin- 
single walker quantum walk. 

Our motivation to use entangled coins in quantum 
walks comes from two sources. First, using entangled 
coins \c) E H n in quantum walks on graphs G(V, E) with 
deg(wi) = m Vi>i € V in which n > m, motivates the 
employment of different shift operators and therefore ex- 
pands the dynamics of the quantum walk. In particular, 
in this paper we use maximally entangled coins in quan- 
tum walks on an infinite line along with shift operators 
with "rest sites" , i.e. states that allow the walker to stay 
at the current vertex. Indeed, it is also possible to in- 
troduce pairs of coins in a classical random walk on an 
infinite line in order to expand its dynamics, but that is 
at the expense of varying the amount of correlation be- 
tween the random variables produced with the outcomes 
of corresponding coins. 

Second, an entangled coin comprised of two qubits, 
each residing in 7i 2 , can be viewed as a single coin de- 
fined on Ji , and then appropriately partitioned. In- 
deed the orthonormal basis {|00), |01), |10), |11)} spans 
the space 7i 4 . However, the phenomenon of entanglement 
represents a supercorrelation between possibly spacelike- 
separated subsystems of a total quantum system. It is 
certainly feasible to generate entanglement between two 
qubits and then separate them either for the purposes 
of an experiment in the laboratory, for example to allow 
for individual addressing of each qubit in some quantum 
information processing experiment (such as implement- 
ing the Hadamard operator), or to send them to opposite 
sides of the universe. This pre-existing entanglement re- 
source is created during a finite time period of interaction 
and then the distinct subsystems can be separated to ar- 
bitrary locations. However, a single four-level quantum 
system cannot be physically broken and spatially sepa- 
rated into two pieces so that each piece is subsequently 
subjected to local operations. Partitioning a single coin 
into two entangled subsystems is nevertheless equivalent 
to using two distinct entangled coins, provided that the 
subsystems do not require to be physically separated for 
practical purposes. Although the mathematical descrip- 
tion is the same, we choose to work with a pair of entan- 
gled qubits. Generation of photonic entangled states, for 
example by way of spontaneous parametric downconver- 
sion, or entangled states in ion traps, is already experi- 
mentally achievable. As such, identifying the entangled 
coins as bipartite states, rather than single an appropri- 
ately partitioned single system residing in H , is a nat- 
ural choice to highlight and motivate possible links to 
experiment. 

Two specific physical implementations of quantum 
walks, namely cavity-QED based ^2 an d ion-trap based 
[43| motivate the scenario considered by us. In both 
these implementations, two-level atoms serve as coins, 
while a cavity mode [42l j or a vibrational mode |43| serve 



as the walker. Two atomic qubits in an ion trap are 
already feasible, and have been prepared in Bell states 
0| and there are several proposals for entangling two 
atomic qubits in a cavity, such as |45| . These atoms can 
then be used as entangled coins with a common cavity 
mode or the common ion trap vibrational mode acting 
as the walker controlled by both these coins. While it 
is straightforward to treat the two atoms as individual 
qubits during this process, it is rather difficult to do en- 
tangling operations between them during the walk with- 
out using/affecting the cavity or vibrational mode which 
is already acting as the walker. Of course, a cavity mode 
or vibrational mode other than the one acting as walker 
may be used to do entangling gates between the atoms, 
but this is complicated. Moreover, in some cases, no 
method of accomplishing a direct unitary entangling gate 
between two atoms may be present, and their initial en- 
tanglement (needed for the walk considered here) may 
have been produced usin g ot her mechanisms (such as de- 
cays and measurements |43). Because of this inherent 
difficulty of doing an entangling gate between the coins 
during a quantum walk, it is easier to imagine a scenario 
of two entangled coin qubits rather than a single four di- 
mensional coin. Once two atoms have been trapped and 
entangled in a cavity, and this has already been done for 
ion traps, the implementation of our scenario is no more 
complex than a single coin quantum walk as the same 
global fields can be applied to both atoms for the coin 
and the shift operations (there is no need for addressing 
the atoms separately). 

In this paper we shall discuss the behaviour of a quan- 
tum walk on an infinite line (also called unrestricted 
quantum walk) with one coin composed of two maximally 
entangled particles, and one walker. We compare the per- 
formance of such a walk with that of a classical random 
walk with one walker and two maximally correlated coins. 
We also show that the use of different shift operators on 
2 and 3 qubit coins leads to different position probability 
distributions in both one and two-dimensional graphs. 

The idea behind correlated coins is simple. For a pair 
of correlated coins C\ and Ci with corresponding out- 
comes [H\,T\) and (fi^,^) one expects that, after ob- 
taining a certain outcome for coin C\, coin Ci will pro- 
duce its corresponding outcome according to a probability 
distribution defined by the degree of correlation between 
both coins. 

For example, the behaviour of a maximally correlated 
pair of coins would be the following: outcomes for coin 
C\ would be given according to a certain probability dis- 
tribution. Let us suppose that coin C\ is unbiased, thus 
outcomes H\ and T\ may each occur with equal probabil- 
ity. Now let us suppose that we get Hi (Ti) as outcome. 
Since the coin pair is maximally correlated, then the out- 
come for coin C2 will certainly be H2 (T2). 

If the degree of correlation were less than maximal be- 
tween coins C\ and C2, then obtaining outcome H\ for 
C\ would imply that the probability of getting H2 as out- 
come for coin C'2 would not be unity. In fact the proba- 
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bility would scale as a monotonically increasing function 
of the degree of correlation between the coins. 

Using correlated coins in classical random walks is 
straightforward. For a classical random walk with a max- 
imally correlated pair of coins it is natural to assign the 
walker one step to the right whenever the pair (Hi, H 2 ) 
(say) is the resulting outcome, and one step to the left 
for the outcome (Ti,T 2 ) (say). In this case, outcomes 
(i?i,T 2 ) and (Ti,H 2 ) have probability zero. 

Indeed, we could enrich our classical random walk by 
allowing coin outcomes 3 = (Hi,T 2 ) and (9 4 = (T 1 ,H 2 ). 
For example, one could use outcome O3 to permit the 
walker to remain in its current position or, alternatively, 
all four outcomes could be used to perform a random walk 
with 1 and 2 steps to the right and left, respectively. In 
particular, the introduction of outcomes that allow rest 
states is a feature used to remove the parity property 
of classical random walks, which consists of finding the 
walker only in even (odd) positions in an even (odd) time 
step. However, the introduction of outcomes O3 and O4 
implies that the coin pair would no longer be maximally 
correlated. 

In the following section we formally introduce a clas- 
sical random walk with a maximally correlated pair of 
coins. In Section III we present a brief background intro- 
duction to unrestricted quantum walks on a line with a 
single coin, followed by our results on unrestricted quan- 
tum walks on a line with a maximally entangled coin. 

Our results show that probability distributions of 
quantum walks with maximally entangled coins have par- 
ticular shapes that are highly invariant to changes in coin 
operators. These results are then compared with those 
obtained for classical random walks with maximally cor- 
related coins. 



II. CLASSICAL RANDOM WALK WITH 2 
MAXIMALLY CORRELATED COINS 



where Cov(X,Y) and Var(X) are the covariance and the 
variance of the corresponding random variables. The 
function p is bounded by — 1 < p < 1. p{C\,C 2 ) = 
means that random variables C\ and C 2 are totally un- 
corrected (i.e. C\ and C 2 are independent), whereas 
p(Ci,C 2 ) = 1 means that random variables C\ and C 2 
are maximally correlated. The case p(Ci,C 2 ) = —1 cor- 
responds to perfect anticorrelation. 

Now consider a classical random walk that has a max- 
imally correlated pair of coins, i.e. p(Ci,C 2 ) — 1- Also 
suppose that the first coin C\ is unbiased. Then, as 
explained in the previous section, the only two out- 
comes allowed for this coin pair are 0\ — {Hi,!!*}) or 
O2 = (Ti, T 2 )- If 0\ allows the walker to move one step 
to the left and 2 allows the walker to move one step 
to the right, it is then clear that using such a coin pair 
in a classical random walk would produce a probability 
distribution equal to that of Eq. ([TJ), with p = h. A plot 
of Eq. JTJ) with number of steps n = 100 and p = i is 
provided in Fig. ([TJ) for the purpose of comparison with 
results presented in Section III. 




A classical result from stochastic processes states that, 
for an unrestricted classical random walk starting at po- 
sition zq = 0, the probability of finding the walker at po- 
sition k after n steps, when with probability p the walker 
takes a step to the right and with probability q = 1 — p 
takes a step to the left (i.e. tossing the coin with prob- 
ability p of obtaining outcome T and probability q of 
obtaining outcome H), is given by 



P 



(n) _ 



±(fc 



k) 



(1) 



for i(fc + n) £ {0,1, ... ,n} and otherwise. 

Tossing a pair of coins produces two discrete random 
variables C\ and C 2 , and the correlation p between these 
two random variables is given by f|36jp 



p(Ci,C 2 ) = 



Cov(Ci,C 2 ) 
^Var(Ci)Var(C 2 ) 



(2) 



FIG. 1: Plot of P, 



(n) 



(i ( ;+>^ +r V ( "- fe) for n 



100 



and p = \ . The probability of finding the walker in position 
k = is equal to 0.0795. Only probabilities corresponding to 
even positions are shown, as odd positions have probability 
equal to zero. 



As can be seen in Fig.Q, the use of maximally cor- 
related unbiased coins in classical random walks is not 
different with respect to a classical random walk with 
a single unbiased coin, as the probability distributions 
from both kinds of classical random walks are exactly 
the same. 

In the following sections we shall compare the results 
obtained by the computation of classical random walks 
with maximally correlated (classical) coins with those of 
quantum walks with maximally entangled (quantum cor- 
related) coins. 
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III. QUANTUM WALKS WITH ENTANGLED 
COINS 

A. Review of Quantum Walks on an Infinite Line 

We now review the mathematical structure of a quan- 
tum walk on a line with one coin and a single walker, with 
a view to using this structure to construct a model for 
unrestricted quantum walks on the line with a maximally 
entangled coin. 

The main components of a quantum walk on a line are 
a walker, a coin, evolution operators for both walker and 
coin, and a set of observables: 

Walker and Coin: The walker is a quantum system liv- 
ing in a Hilbert space of infinite but countable dimension 
Hp. It is customary to use vectors from the canonical 
(computational) basis of H p as "position sites" for the 
walker. So, we denote the walker as |position) G H p and 
affirm that the canonical basis states \i) p that span H p , as 
well as any superposition of the form ^\ a iV)p subject to 
Tlii \ a i\ 2 — 1> are valid states for | position). The walker 
is usually initialized at the 'origin', i.e. | position) initial = 
|0) P . 

The coin is a quantum system living in a 2-dimensional 
Hilbert space H c . The coin may take the canonical ba- 
sis states |0) and |1) as well as any superposition of 
these basis states. Therefore |coin) e H c and a general 
normalized state of the coin may be written as |coin) 
= a\0) c + b\l) c , where \a\ 2 + \b\ 2 = 1. 

The total state of the quantum walk resides in Ht — 
H p <8> H c - Only product states of Ht have been used 
as initial states, that is, ^initial = | position) initial <8> 

| COin) initial- 

Evolution Operators: The evolution of a quantum 
walk is divided into two parts that closely resemble the 
behaviour of a classical random walk. In the classical 
case, chance plays a key role in the evolution of the sys- 
tem. This is evident in the following example: we first 
toss a coin (either biased or unbiased) and then, depend- 
ing on the coin outcome, the walker moves one step either 
to the right or to the left. 

In the quantum case, the equivalent of the previous 
process is to apply an evolution operator to the coin state 
followed by a conditional shift operator to the total sys- 
tem. The purpose of the coin operator is to render the 
coin state in a superposition, and the randomness is in- 
troduced by performing a measurement on the system 
after both evolution operators have been applied to the 
total quantum system several times. 

Among coin operators, customarily denoted by C, the 
Hadamard operator 

H = -J=(|0) cc (0| + |0) cc (l| + |l) cc (0| - |l) cc (l|) (3) 

has been extensively used. 

For the conditional shift operator use is made of a uni- 
tary operator that allows the walker to go one step for- 
ward if the accompanying coin state is one of the two 



basis states (e.g. |0)), or one step backwards if the ac- 
companying coin state is the other basis state (|1)). A 
suitable conditional shift operator has the form 

S= |0) cc <0|®^|i + l) pp <i| 

i 

+ |l}cc(l|®]>>- 1 U J l- ( 4 ) 

i 

Consequently, the operator on the total Hilbert 
space is U = S.(C ® l p ) and a succint mathemati- 
cal representation of a quantum walk after n steps is 

|V>) = (17) n |l/>)initial, Where ^initial = I position) initial ® 
I COin) initial- 

Observables: The advantages of quantum walks over 
classical random walks are a consequence of interference 
effects between coin and walker after several applications 
of U . However, we must perform a measurement at some 
point in order to know the outcome of our walk. To do 
so, we define a set of observables according to the basis 
states that have been used to define coin and walker. 

In order to extract information from the composite 
quantum system, we first perform a measurement on the 
coin using the observable 

M c = a |0)cc(0|+ai|l)cc(l|. (5) 

A measurement must then be performed on the position 
states of the walker by using the operator 

M p = ^2ai\i) pp (i\. (6) 



B. Mathematical Structure of Quantum Walks on 
an Infinite Line Using a Maximally Entangled Coin 

As before, the elements of an unrestricted quantum 
walk on a line are a walker, a coin, evolution operators for 
both coin and walker, and a set of observables. We shall 
provide a detailed description of each element motivated 
by the previous subsection. 

Walker and Coin: The walker is, as in the unre- 
stricted quantum walk with a single coin, a quantum 
system | position) residing in a Hilbert space of infinite 
but countable dimension Hp. The canonical basis states 
\i)p that span Hp, as well as any superposition of the 
form Y^i a i\i)p subject to J2i \ a i\ 2 = 1> are valid states 
for the walker. The walker is usually initialized at the 
'origin' i.e. |position) = \0)p. 

The coin is now an entangled system of two qubits i.e. 
a quantum system living in a 4-dimensional Hilbert space 
Hec- We denote coin initial states as |coin) - Also, we 
shall use the following Bell states as coin initial states 

|*+) = -±=(|00) + |11)) (7a) 
|<&-) = 4(100} (7b) 
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|*+) = _(|01) + |10)) 



(7c) 



which are maximally entangled pure bipartite states with 
reduced von Neumann entropy equal to unity. We shall 
examine the consequences of employing such maximally 
entangled states by comparing the resulting walks with 
those resulting from using maximally correlated coins in 
classical random walks. The Bell singlet state \^^) = 
^7g(|01) — 1 1 0) ) is not employed as an entangled coin as 
it is left invariant when the same local unitary operator 
is applied to both coins. 

The total initial state of the quantum walk resides in 
the Hubert space Ht = Hp <£> TLec and has the form 

|V>)o = |position) ® |coin) (8) 

Entanglement measure: In order to quantify the de- 
gree of entanglement of the coins used in this paper, we 
shall employ the reduced von Neumann entropy mea- 
sure. For a pure quantum state of a composite 
system AB with dim(A) = d,A and dim(B) = ds, let 
IV') = Z)i=i a »lu>|iB>, (d = mm(d A ,d B ), a t > and 
Sj=i a i = 1) be its Schmidt decomposition. Also, let 
Pa = trs(|V')(V , |) and Pb = tr^dV") (V 1 !) be the reduced 
density operators of systems A and B respectively. The 
entropy of entanglement E(\ip)) is the von Neumann en- 
tropy of the reduced density operator [33] 



E(m = S{ PA ) = S( PB ) 



i=l 



a?log 2 (a?). (9) 



E is a monotonically-increasing function of the entangle- 
ment present in the system AB. A non-entangled state 
has E = 0. States \ip) e H d for which E(ip) = d are called 
maximally entangled states in d dimensions. In particu- 
lar, note that for those quantum states described by Eqs. 

0, JTbJ and O E(\<P+)) = J5(|* - )) = E(\f+)) = 1, 

1. e. these states are maximally entangled. 
Evolution Operators: The evolution operators used 
are more complex than those for quantum walks with 
single coins. As in the single coin case, the only require- 
ment evolution operators must fulfil is that of unitarity. 

Let us start by defining evolution operators for an en- 
tangled coin. Since the coin is a bipartite system, its 
evolution operator is defined as the tensor product of 
two single-qubit coin operators: 

C EC = C®C (10) 

For example, we could define the operator C§ c as the 
tensor product H® 2 : 



|01)(00| + |10)(00| + |11)(00| 

+ |00)(01|-|01){01| + |10)(01|-|11)(01| 
+ |00)(10| + |01)(10|-|10)(10|-|11)(10| 
+ |00)(11|-|01)<11|-|10)(11| + |11)(11|). 



An alternative bipartite coin operator is produced 
by computing the tensor product Y® 2 where Y — 
^(|0)(Q|+<|0)(1|+*|1)(0| + |1)<1|), namely 



C Y 



EC 



:-(|00)(00M 

+ t|00)(01| 
+ t|00)(10| 
-|00)(11|- 



z|01>(00| ■ 

f |01)(01| 
-|01)(10| 
■i|01)<ll| 



i|10)(00| 

•|10)(01|- 
-|10)(10|- 
-t|10)<ll| 



-|11><00| 

i|ll>(01| 
i|ll>(10| 

HH)(U|). 

(12) 



Both coin operators are fully separable, thus any en- 
tanglement in the coins is due to the initial states used. 

The conditional shift operator Sec necessarily allows 
the walker to move either forwards or backwards along 
the line, depending on the state of the coin. The operator 



Sec = |00) cc (00|®J]|z + l) 

i 

+|01} co (01|®$3|i> 

i 

+|10) ec (10|®5^| 

i 

+in> cc (iii®5^ii-i) H ,(i| 



ipp\ 



(13) 



embodies the stochastic behaviour of a classical random 
walk with a maximally correlated coin pair. It is only 
when both coins reside in the |00) or |11) state that the 
walker moves either forwards or backwards along the line; 
otherwise the walker does not move. 

Note that Sec is one of a family of valid definable shift 
operators. Indeed, it might be troublesome to identify a 
classical counterpart for some of these operators: their 
existence is uniquely quantum-mechanical in origin. One 
such alternative operator is 



S'ec = |00) cc (00|®^|z + 2) 

i 

+|01) cc (01|®J]|i+l 

i 

+|io) cc <io|®^K-i 

i 

+|ll) ec <ll|(8>5^|i-2) 



vi>\ 



vi>\ 



(14) 



(11) 



VP \ 



The total evolution operator has the structure Ut — 
Sec-{Cec ® Ip) and a succint mathematical representa- 
tion of a quantum walk after N steps is = (Ut) n \iP)o> 
where \ip)o denotes the initial state of the walker and the 
coin. 

Observables: The observables defined here are used to 
extract information about the state of the quantum walk 
|V) = (Ut) n \^)o- 
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We first perform measurements on the coin using the 
observable 

M EC = A)o|00) cc (00| + /3oi|01) cc (01| 

0io|io)cc<io| + Ai|n>cc<n|. 

Measurements are then performed on the position states 
using the operator 

Mr ( 16 ) 

j 

With the purpose of introducing the results presented 
in the rest of this paper we compare in Table 1 the actual 
position probability values for a classical random walk on 
an infinite line (Eq. Q), and a quantum walk with initial 
state |<i> + ) = -7= (| 00) + 111)) and coin and shift operators 

given by Eqs. and (JT3J, respectively. 

Table 1. Position Probability values for classical 
random walk and quantum walk 
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FIG. 2: For both plots, coin initial state is |$+) = ^jdOO) + 
1 11}) and the number of steps is 100. Coin operators for red 
and dotted blue plots are given by Eq. JTTJ and Eq. itHl 
respectively. 
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C. Results for Quantum Walks on an Infinite Line 
Using a Maximally Entangled Coin 
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FIG. 3: For both plots, coin initial state is |$ + } = ^jdOO) + 
1 11}) and the number of steps is 200. Coin operators for red 
and dotted blue plots are given by Eq. (JTTJ and Eq. ifHSjl 
respectively. 



In order to investigate the properties of unrestricted 
quantum walks with entangled coins, we have computed 
several simulations using bipartite maximally entangled 
coin states described by Eqs. J7aJ, (|7b|) and lj7c|) . and 
coin evolution operators described by Eqs. ltTT|) and ((T2^l . 
In all cases, the initial position state of the walker cor- 
responds to the origin, i.e. |position) = |0). The shift 
operator employed is, with the exception of Fig. JSJ, that 
of Eq. (|T3J. 

Let us first discuss the quantum walks whose graphs 
are shown in Fig. The initial entangled coin state is 
given by Eq. (|7a|l and the number of steps is 100. For the 
red plot in Fig. @ the coin operator is given by Eq. (|llfl . 
while for the dotted blue plot in the same Fig. @ the 
coin operator is that of Eq. I|12|) ■ 

The first notable property of these quantum walks is 
that, unlike the classical case in which the most probable 
location of the walker is at the origin and the probability 



distribution has a single peak, in the quantum case a cer- 
tain range of very likely positions about the position |0) 
is evident but in addition there are a further two regions 
at the extreme zones of the walk in which it is likely to 
find the particle. This is the 'three peak zones' property 
of the shift operator defined in this way. The 'three peak 
zones' property could mean an additional advantage of 
quantum walks over classical random walks. 

We also note that the probability of finding the walker 
in the most likely position, |0), is much higher in the 
quantum case (~ 0.171242 in red plot of Fig. (0 and 
~ 0.221622 in dotted blue plot of Fig. J3J)) than in the 
classical case (~ 0.0795). Incidentally, we find that the 
use of different coin initial states maintains the basic 
structure of the probability distribution, unlike the quan- 
tum walk with a single coin in which the use of different 
coin initial states can lead to different probability distri- 
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butions ([3, [13 and H3) . 

The position probability distributions shown in Fig. iffil 
could embody some advantages when used in an appro- 
priate application framework. For example, let us sup- 
pose we want to design algorithms whose purpose is to 
find the wrong values in a proposed solution of a prob- 
lem (a concrete case is to find all wrong binary values 
assigned to the initial conditions of the algorithm used 
in [38j to find a solution to the 3-SAT problem.) We use a 
100-steps classical random walk (Fig. QJ) to design algo- 
rithm C and a 100-steps quantum walk with maximally 
entangled coins (red plot of Fig. iffil) to design algorithm 
Q. 

Depending on the actual number of wrong values, the 
probability distribution of red plot of Fig. J5J) could 
help to make algorithm Q faster than algorithm C. For 
example, if the number of wrong values is in the range 
40 - 70, the probability of finding the quantum walker 
of Fig. (|2J) is much higher than finding the classical 
walker of Fig. Q. Actual probability values (note the 
differences in orders of magnitude) are shown in Table 
2. It must be noted that employing a quantum walk on 
a line with a single coin for building algorithm Q would 
also produce higher probability values than a classical 
random walk in those positions shown in Table 2, thus 
the choice of quantum walk could depend on some other 
factors like implementation feasibility. 

Table 2. Position Probabilities for Classical and 
Quantum Walkers 



Position 


Classical Walker 


Quantum Walker 


40 


2.31 x 10~ 5 


1.80 x 10~ 3 


50 


1.91 x 10~ 7 


1.50 x 10~ 3 


60 


4.22 x lO" 10 


1.03 x 10~ 2 


70 


1.99 x 10~ 13 


3.78 x 10~ 2 



-200 -150 -100 



100 150 200 



FIG. 4: Coin initial state is | $") = -^(|00)-|11». Number of 
steps is 100 for red and dotted blue plots and 200 for starred 
yellow and black dotted plots. Coin operator for red and 
starred yellow plots is given by Eq. I)ll|l . while coin operators 
for dotted blue and dotted black plots is given by Eq. (II 2t 
respectively. 



0.16 
0.14 
0.12 
0.1 
0.08 
, 0.06 



J, 



-50 50 

Position 



A consequence of the previous two properties of the 
quantum walk is a sharper and narrower peak in the 
probability distribution around position |0). Again, this 
may be of some advantage depending on the application 
of the quantum walk (for example, less dispersion around 
the most likely solution to the computational problem 
posed in the two previous paragraphs). 

The probability distributions for quantum walks in 
Fig. are very similar in structure to those of Fig. J5J, 
the only difference being the number of steps (200 as op- 
posed to 100). For Fig. J2J the initial entangled coin state 
is given by Eq. Q7af) . Eq. Qllf) is used as the coin operator 
in the red plot of Fig. (J3J) whereas Eq. (|12fl is the coin 
operator for the dotted blue plot of Fig. J2J). For 200 
steps the peaks on both extreme zones are smaller than 
for 100 steps, the reason being the increased number of 
small probabilities that correspond to those regions be- 
tween the extreme peaks and the central peak. A wider 
region is covered in the case of 200 steps than for 100 
steps. 



FIG. 5: Coin initial state is |*+) = -i (|01) + 110>). Number of 
steps is 100 for red and dotted blue plots and 200 for starred 
yellow and black dotted plots. Coin operator for red and 
starred yellow plots is given by Eq. I|ll|l . while coin operators 
for dotted blue and dotted black plots is given by Eq. W'll 
respectively. 



An examination of Figs. @ and © is straightfor- 
ward, as their bulk properties closely resembling those 
of Fig. J2J and Fig. JSJ. The probability distributions in 
Fig. and were computed using Eq. (|7b|l as the initial 
coin state and the same initial conditions and shift oper- 
ators as for Fig. © and Fig. ©. The number of steps 
is 100 for red and dotted blue plots and 200 for starred 
yellow and black dotted plots. Coin operator for red and 
starred yellow plots is given by Eq. (|llf> . while coin op- 
erators for dotted blue and dotted black plots is given by 
Eq. (|12J) respectively. The 'three-peak zones' feature is 
again evident. Furthermore, the bulk properties of the 
probability distributions are highly invariant to changes 
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Entangled quantum walk vs. 2-step one coin quantum walk 
0.08 | 1 1 1 1 1 1 ,— 

0.07 - 

0.06 : ; :• 

0.05 
1 0.04 - 

s 

0.03 - 




-200 -150 -100 -50 50 100 150 200 

Position 

FIG. 6: For the red line plot, the coin initial state is given 
by |* + ) = TTjdOO} + |H})- Coin operator is given by Eq. 
llll l and shift operator by Eq. H14C . For the blue dashed 
graph, coin initial state is given by (V0.85[0) c — V0.15|l) c ), 
coin operator is the Hadamard operator (Eq. 0) as coin 
operator and shift operator given by |0}(0| ® Yli I* + 2)(i| + 
® Si N — In both cases, the number of steps is 

100. 

in coin operators (there is a slight difference in the prob- 
ability distribution value at the origin). In both cases 
the probability distribution value at the origin is much 
larger than in the classical random walk case. A similar 
discussion having Eq. JTcJ) as coin initial state, and using 
the same colors as in Fig. (@J to refer to coin operators 
and number of steps, applies to Fig. (JSJ. 

In order to further motivate the richness of quantum 
walks with entangled coins, we present the graph shown 
in Fig. 10 red line plot) computed using Eq. l(7a|) as the 
initial state of the coin, Eq. (|llfl as the coin operator and 
Eq. I|14|) as the shift operator. This graph closely resem- 
bles that of a 2-step quantum walk Fig. JfjJ dotted blue 
plot) with initial state (V0M\0) c - vCl5|l> c )® |0> p pj . 
Hadamard operator (Eq. (J3J) as coin operator and shift 
operator given by |Q)<0|®£)< |*+2)(*|+|l)(l|®Ei K~ 2> <i| 
(the number of steps in both walks is 100). However, the 
graph corresponding to the quantum walks with a maxi- 
mally entangled coin has no parity restriction, as opposed 
to the 2-step quantum walk, and this explains the higher 
probability values for the 2-step quantum walks. 

As opposed to the previous cases (Figs. 121-03) in which 
the walker remains static when the quantum coin state 
component is either |01) or |10), in this case the walker is 
forced to jump either one or two steps, depending on the 
components of the coin state. As can be seen in Fig. JBJ, 
the behaviour of the quantum walk dramatically changes 
as a consequence of the change in the shift operator. In 
this case, constructive interference takes place not only in 
certain areas of the graph (as is the case with the 'three 
peak zones' property) but in a wider region. Indeed, this 
walk bears a resemblance to a quantum walk using a 
single walker and a single Hadamard coin p2^ . 

Finally we would like to emphasize that in stark con- 



trast to the probability distributions of the classical case 
in which only certain walker positions have a probability 
different from zero, namely those positions whose parity 
is that of the total number of steps, in the quantum cases 
presented in this paper we observe no such constraint on 
the numerical data produced. As stated in the introduc- 
tion, the dynamics of classical random walks can remove 
the parity constraint by permitting the use of 'rest sites' 
at the expense of varying the amount of correlation be- 
tween the coins. 



IV. QUANTUM WALKS USING COINS WITH 
DIFFERENT VALUES OF ENTANGLEMENT 

In order to compare the properties of quantum walks 
with coins having different degrees of entanglement, we 
present in this section several probability distributions 
computed using bipartite coins. The graphs of those 
probability distributions are shown in Figs. J3-|3J). All 
graphs shown in Figs. {7|-|H1) were computed using Eq. 

as coin operator and Eq. H13[) as shift operator. The 
initial position state in all cases is the origin, i.e. |0) p . 

The probability distribution presented in Fig. {7J) 
shows a typical skewed (asymmetric) behaviour in quan- 
tum walks. The graph is produced using the bipartite 
quantum state \9q) = ^(|0) + |1))(|0) + |1)) as coin initial 
state (fromEq. JJjJ, E(\9 )) — so \9 ) is non-entangled). 
This graph resembles the behaviour of a quantum walk 
presented in |29j using a coin in initial state |00) {\RR) 
in their notation). 

Let us now focus on the behaviour of the quantum walk 
shown in Fig. (|SJ, which was produced using a partially 
entangled initial coin state. The coin was initialized in 

the state = ±|00) + ±|01) + ^±|10) + ^±±|11). 
Again, using Eq. ©, we find that E(\6i)) = 0.5, i.e. |6>i) 
is partially entangled. 

We can see that an immediate effect of an entangled 
coin initial state is the development of a third peak, in 
the case a peak on the LHS of the graph. This third 
peak reduces the skewness of the probability distribution 
computed with a non-entangled coin initial state (Fig. 
©)■ 

Let us now compare Figs. (0 and (|10f) with the red 
plot from Fig. @, created with the maximally entan- 
gled state -j)j(|00) + 1 1 1 ) ) as initial coin state. We can 
see that the increasing use of entanglement provides a 
greater degree of symmetry to the resulting probability 
distribution. 

Now, if we expand the properties of initial conditions 
by allowing coins to be initialized in states with complex 
coefficients, we would then obtain probability distribu- 
tions that would be similar to those of quantum walks 
with maximally entangled coins with real coefficients. 

For example, we show in Fig. © the position probabil- 
ity distribution computed with initial coin state |(| 0) + 
i|l))(|0) + *|1))- Fig. @ bears a striking resemblance 
to those of Figs. (0-EJ). However, even though qualita- 
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FIG. 7: Quantum walk computed with a coin initialized in the 
state |(|0) + |1))(|0) + |1)), i.e. a non-entangled state with 
real coefficients. 100 steps, coin and shift operators given by 
Eqs. 1111 and 1131 respectively. 



FIG. 8: Quantum walk computed with a coin initialized in the 
state || 00) + §|01) + ^%±|10) + ^±i|ll), i.e. a partially- 
entangled state with real coefficients. 100 steps, coin and 
shift operators given by Eqs. il l H and l|13|l respectively. Note 
that the entanglement of the coin initial state reduces the 
assymetry of the graph by creating a new third peak. 



tively a three peaked structure is evident, quantitatively 
it differs considerably. To illustrate this numerical dif- 
ference, we appeal to Fig. (JTDJ, which compares three 
different quantum walks. The probability distribution 
computed with coin -i=(|00) + |11)) corresponds to the 
blue starred points on the graph, while the probability 
distributions computed with coins -i=(|00) — |11)) and 

1 1 00) + ||01) + || 10) - ||11) are depicted using red dots 
and yellow circles respectively. Thus the entanglement 
of the initial coin state helps to both tune up and tune 
down the ratio of the central peak to the side peaks. 

Numerical values show that entanglement plays an ac- 
tive role in the actual probability of finding the walker 
in a certain position. For example, consider walker po- 



FIG. 9: Coin initial state is |(|0) +i|l))(|0) + i|l)). 100 steps, 
coin and shift operators given by Eqs. il l U and il'At respec- 
tively. The use of complex coefficients in the coin initial state 
delivers a symmetric probability distribution very similar to 
those shown in Figs. 121-151. 




FIG. 10: The probability distribution computed with coin 
■^75 (1 00) + 1 11)) corresponds to the blue starred graph. Prob- 
ability distributions computed with coins -^=(|00) — |11)) and 

||00) + ||01) + 1 1 10) - 1 1 11) are shown in red dots and yel- 
low circles respectively. All graphs were computed after 100 
steps using Eq. HI IB as coin operator and Eq. 1131 as shift 
operator. 



sitions 60-70. The highest values in this region are at- 
tained by the probability distribution computed with 
coin 771 (|00) + 1 11)). In a different region, that of the 
central peak, the probability distribution of the non- 
entangled coin initial state ||00) + ||01) + ||10) - ||11) 
is between those values produced by the two probabil- 
ity distributions obtained by computing quantum walks 
with maximally entangled states. 

Another example of a symmetric graph produced using 



coins in non-entangled states with complex coefficients 
has been presented by Inui and Konno in |46j . This sym- 
metric graph was produced using a coin initialized in the 
state 1 1 00) + ||01) + ||10> + ||11>. 



V. QUANTUM WALKS WITH MORE THAN 
TWO MAXIMALLY ENTANGLED COINS 



An interesting property of using several entangled 
qubits as coins is the fact that the number of coin and 
shift operators available for use also increases. Conse- 
quently, several different position probability distribu- 
tions can be computed. 

For example, in Fig. Hill red plot) the graph of a 100- 
steps quantum walk with the GHZ state -^=(|000) + |111)) 
as coin initial state is shown, the coin operator being 
given by H® 3 where H is Hadamard operator (Eq. 
and shift operator given by 



S 3a = |000) cc (000| ® N + Vppi* 

i 

+ |001) cc (001|®^|z) pp (« 

i 

+|Q10)cc<010| ® !*)»<* 

i 

+ |011) cc (011|<8>^|i>p P (i 

i 

+ |1OO) cc (1OO|0^|z) pp (« 

i 

+|ioi)cc<ioi|®£) 

i 

+ |110) cc (110|<8>^|i>p P (i 

i 

+|m) cc (m|®^|i- i) pp (i 



which has a 4-peak probability distribution. However, 
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changing the shift operator to 

S 3b = |000) cc (000|®^|i + 3) pp (z| 

i 

+ |001) cc (001|®J]|i + 2) pp (z| 

i 

+1010)00(0101 ® 1* + VppW 

i 

+|011) cc <011|®5^|i) w <»| 

(18) 

+ |100) cc (100|®^|z) pp (z| 

i 

+|101) CC (101|®5Z|*-1) PP (<| 

i 

+|110) cc (110|®^|i-2) w (i| 

i 

+|lll) cc (lll|®5^|i-3) w <i| 

i 

results in the blue plot of Fig. (|ll|l which has no such 
readily evident peak structure. 

0.06 1 1 1 1 1 1 1 1 

I 0.05- 
o 



Q. 

To 0.04- 




300 



Position 

FIG. 11: Position probability distributions for two quantum 
walks on a line with GHZ state -^=(|000) + |111)) as initial 

tripartite coin state and coin operator //® 3 . The red plot 
was computed using the shift operator in Eq. 1171 and the 
blue plot using the shift operator in Eq. 1181 . While the red 
plot shows an evident 4-peak structure, the blue plot does not 
present such a behaviour. 

The potential richness of quantum walks increases 
when taking into consideration graphs of more than one 
dimension (efforts to understand the properties of quan- 
tum walks on graphs are presented in |l3j . while a pro- 
posal for a physical realization of a 2-dimensional quan- 
tum walk is given in jUj). For example, Fig. Ijl2(l shows 
the peak structure of a 50-step quantum walk on a graph 
with initial state given again by 4=(|000) + 1 1 1 1) ) , coin 
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operator given by H® 3 and shift operator given by Eq. 



Sec = |000)cc<000|®5^|* + l > j) w <*,j| 

i 

+ 1001) i j)pp(h j\ 

i 

i 

+|011) cc (011| ® Yl I* 

(19) 

+ |100) cc (100|®^|z 1 j)pp(h j\ 
i 

+iioi> eo <ioi|®£;ii,i-i) 

ppV") 3\ 

i 

+|110) cc (110|®5^|i > i) sp <<,j| 

i 

+iiii) cc (ini®5^i*-i,i) 



Quantum walk on a graph 
with coin |c> = (1/2 1/2 )(|000>+|111» 



0.012^ 
0.01 - 
& 0.008 - 




i Position SU " 50 i position 



FIG. 12: Position probability distribution of a quantum walk 
on a 2-dimensional graph computed with coin initial state 
^=(|000) + | 111)) and shift operator given by Eq. JTHJ. The 
number of steps is 50. The graph has 2 high peak regions and 
several other small peaks in the central region. 
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